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Q . ABSTRACT 

O 

The current status of extragalactic deuterium abundance is discussed using two examples 
of 'low' and 'high' D/H measurements. We show that the discordance of these two types of 

^ ■ D abundances may be a consequence of the spatial correlations in the stochastic velocity field. 

cn . Within the framework of the generalized procedure (accounting for such effects) one finds good 

agreement between different observations and the theoretical predictions for standard big bang 
nucleosynthesis (SBBN). In particular, we show that the deuterium absorption seen at Za = 2.504 

r-- ; toward Q 1009+2956 and the H+D Lya profile observed at Za = 0.701 toward Q 1718+4807 are 

5^ I compatible with D/H -^ 4.1 — 4.6 x 10^^. This result supports SBBN and, thus, no inhomogeneity 

1-^ ■ is needed. The problem of precise D/H measurements is discussed. 

O ' Introduction 

j^ I The measurement of deuterium abundance at high redshift from absorption spectra of QSOs 

is the most sensitive test of physical conditions in the early universe just after the era of big 
bang nucleosynthesis. The standard BBN model predicts strong dependence of the primordial 
ratio of deuterium to hydrogen nuclei D/H on the cosmological baryon-to-photon ratio ry and 

C^ ' the effective number of light neutrino species Ni, (e.g. Sarkar 1996). According to the basic idea 

of homogeneity and isotropy of big bang theory the primordial deuterium abundance should 
not vary in space. One can only expect that the D/H ratio decreases with cosmic time due to 
conversion of D into ^He and heavier elements in stars. 

On the other hand, current measurements of D abundance at high redshift reveal a scatter 
of the D/H values over approximately one order of magnitude. For instance, recent HST obser- 
vations of the Za = 0.701 absorption-line system toward the quasar Q 1718+4807 (Webb et al, 
1997a, b) show D/H =1.8 — 3.1 x 10~^. This hydrogen isotopic ratio is significantly higher than 
that derived from other quasar spectra at Za = 2.504 [D/H =1.8 — 3.5 x 10~^ by Buries &: Tytler, 
1996; D/H = 2.9 - 4.6 x 10"^ by Levshakov, Kegel & Takahara, 1997 (LKT, hereinafter)], and 
at Za = 3.572 [D/H = 1.7 - 2.9 x 10"^ by Tytler et al, 1996; D/H > 4 x 10"^ by Songaila et 
al, 1997]. 
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The apparent spread of the D/H values makes some authors to assume fluctuations in the 
baryon-to-photon ratio at the epoch of BBN, i.e. to assume that big-bang nucleosynthesis has 
occurred inhomogeneously (see e.g. Webb et al., and references cited therein). 

Current modeling of the evolution of D/H does not allow to estimate the primordial deu- 
terium abundance with a sufficient accuracy. The depletion factor is strongly model dependent 
and changes from about 2-3 for the chemical evolution models assuming inflow (e.g. Prantzos 
1996, Tosi et al. 1997) up to 10 for the models employing a galactic wind, i.e. outflow (e.g. Scully 
et al. 1997). This means that both 'high' (~ 2 x 10""^) and 'low' (~ 3 x 10"^) D/H ratios do 
not contradict the SBBN prediction if the ISM deuterium abundance of 1.6 x 10^^ (e.g. Linsky 
&: Wood, 1996) is adopted. Nevertheless, the precise measurements of absolute values of D/H 
at high redshift are extremely important to check whether BBN does occur inhomogeneously 
or the basic idea of homogeneity and isotropy is still valid. The fundamental character of this 
cosmological test requires an unambiguous interpretation of spectral observations. 

Ho^v accurately can be determined D/H ? 

It is well known that the physical parameters derived from spectral data depend on the 
assumptions made with respect to the line broadening mechanism. For interstellar (intergalactic) 
absorption lines a 'non-thermal broadening' is usually assumed to be caused by a large scale 
motion of the absorbing gas. The commonly used microturbulent approach considers bulk 
motions being completely uncorrelated. This yields a symmetrical (Gaussian) distribution of 
the velocity components parallel to the line of sight, which in turn leads to a symmetrical line 
profile function being a convolution of a Voigt profile and a Gaussian. 

Actually, any turbulent flow exhibits an immanent structure in which the velocities in neigh- 
boring volume elements are correlated with each other. Different aspects of the line formation 
processes in correlated turbulent media have been recently discussed in a series of papers by Lev- 
shakov & Kegel (1997, LK hereinafter), Levshakov, Kegel & Mazets (1997, LKM hereinafter), 
and by LKT. 

Once the statistical properties of a turbulent cloud have been specified, the way spectral 
lines ought to be calculated, depends on the problem considered. Considering emission lines one 
is dealing with many lines of sight and the observed intensity should closely correspond to the 
theoretical expectation value (see e.g. Albrecht & Kegel 1987). If, however, one observes the 
absorption spectrum in the light of a point-like background source, one is dealing with one line 
of sight only, and the actually observed intensity may deviate substantially from the expectation 
value (for details see LK and LKM), since averaging along one line of sight only, corresponds 
to averaging over an incomplete sample. The observed line profile is determined by the velocity 
distribution along the particular line of sight. For large values of the ratio of the cloud thickness 
L to the correlation length I the distribution function p{v) for the velocity component along the 
line of sight approaches the statistical average, which has been assumed to be a Gaussian. For 
values of L/l of only a few, however, p(v) may deviate substantially from a Gaussian, and is 
asymmetric in general. This leads to a complex shape of the absorption coefficient for which 
the assumption of Voigt profiles could be extremely misleading. The actual D/H ratio may turn 
out to be higher or lower than the value obtained from the Voigt-fitting procedure. 



The present paper is primarily aimed at the inverse problem in the analysis of the H+D 
Lya absorption observed by Buries & Tytler (1996) and by Webb et al. ( 1997a, b). The original 
analysis was performed in the framework of the microturbulent model. Here we make an attempt 
to re-analyze the observational data on the basis of a more general mesoturbulent model. We 
consider a cloud with a turbulent velocity field with finite correlation length but of homogeneous 
(Hl-)density and temperature. The velocity field is characterized by its rms amplitude at and 
its correlation length /. The model is identical to that of LKT. - The objective is to investigate 
whether the data in question may be interpreted by a unique D/H ratio consistent with the 
primordial abundances of '^He and ^Li and theoretical SBBN predictions. 

Parameter estimation by the RMC procedure 

To estimate physical parameters and an appropriate velocity field structure along the line of 
sight, we used a Reverse Monte Carlo [RMC] technique, i.e. a stochastic optimization algorithm 
developed to solve optimization problems with a very large number of free parameters (see LKT). 

The algorithm requires to define a simulation box for the 5 physical parameters : N(HI), 
D/H, Tkin, ct/vth, and L/l (here Vfh denotes the thermal width of hydrogen lines). ~ The 
continuous random function of the coordinate v{s) is represented by its sampled values at equal 
space intervals As, i.e. by {fi, f2, . . . ,Vk}, the vector of the velocity components parallel to the 
line of sight at the spatial points Sj. The total number of intervals depends on the values of 
at/vth and L/l, being typically ~ 100 for hydrogen absorption lines. 

The convergence of the RMC procedure depends on the size of the simulation box which is 
unknown a priori. If it is too large, the computing time increases considerably, and if it is too 
small, the required minimum of the objective function may not be reached at all. Therefore the 
box size must be chosen with some care and should be adjusted to the experimental data. In 
the next sections we consider two examples of the RMC application. 

Deuterium abundance ad z = 2.504 toward Q 1009-(-2956 

We applied at first the RMC procedure to a template H+D Lya profile which reproduces the 
Q1009+2956 spectrum with the DI Lya line seen at Za = 2.504 (Buries k Tytler 1996). Here 
we present only a part of our results, the full analysis is given by LKT. 

To simulate real data, we added the experimental uncertainties to the template intensities 
which were sampled in equidistant bins as shown in Fig. 1 by dots and corresponding error 
bars. One-component mesoturbulent model with a homogeneous density and temperature was 
adopted. Adequate profile fits for three different sets of parameters are shown in panels (al, 
bl, cl) by solid curves, whereas the individual HI and DI profiles are the dashed and dotted 
curves, respectively. The estimated parameters (D/H, A^hi, ^femj f^t/vth, L/l) and x^un values 
per degree of freedom are also listed in these panels for each RMC solution. 

The essential difference between the results of the two-component microturbulent model 
adopted by Buries & Tytler and ours lies in the estimation of the hydro dynamical velocities 
in the Za = 2.504 absorbing region. It is generally believed that absorption line systems with 
N(HI)'^ 10^^ cm^^ arise in the halos of putative intervening galaxies (e.g. Tytler et al. 1996). 
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Figure 1: (al, bl, cl) - A template H+D Lya profile (dots with error bars) representing the 
normalized intensities and their uncertainties of the z = 2.504 system toward Q 1009+2956 in 
accord with Buries & Tytler (1996). The solid curves show the results of the RMC minimization, 
whereas the dotted and dashed curves are the separate profiles of DI and HI, respectively. Also 
shown are the best fitting parameters and the x^ values obtained for each case. 
(a2, b2, c2) - The corresponding individual realizations of the velocity distribution u{x) in units 

of CJi. 

(a3, b3, c3) - The histograms are the projected velocity distributions p{u). 



The RMC procedure yields at — 25 km s ^, whereas the too-component model leads to 0"^ ~ 2 
km s~^ which is evidently too low as compared with direct observations of galactic halos at 
z> 2 (van Ojik et al. 1997) which show that ctj ~ 40 ± 15 km s-\ if Tfcj„ ~ 10^ K. 

In this particular absorption system, the study of the H+D Lya profile yields D/H = (3.75 it 
0.85) X 10^^ which is shghtly higher than the value 2.5ltQg9 x 10^^ found by Buries & Tytler. 
This difference is, however, very significant because it leads to limits on D/H consistent with 
other observations discussed below. 

Deuterium abundance at z = 0.701 toward Q 1718+4807 

The absorbing material at z = 0.701 provides an accurate determination of the total hydrogen 
column density due to the extreme sharpness of the Lyman break measured by the International 
Ultraviolet Explorer (lUE) satellite (see Webb et al.). 

To analyze this system we adopted for the physical parameters the following boundaries 
(details are given in Levshakov, Kegel & Takahara 1998) : 
1.70x10^^ cm-2 < N(HI) < 1.78 x 10^^ cm^^ 

For D/H we use the range from 3.0x10^^ to 5.0x10^^, trying to find a low D/H solution. 
For Tfcj„ we use the interval 10^ — 2 x 10^ K, and, thus, (Tt/vth uiay range within 1.3 - 4.3. 
For L/l ^ 1 the meso- and microturbulent profiles tend to be identical (see LK). Different 
microturbulent models have been thoroughly investigated by Webb et al., therefore we will 
consider moderate L/l ratios in the range 1.0 - 5.0. 

Following Webb et al. (1997b), we exclude the Si III line from our analysis of the H+D Lya 
profile. As shown by Vidal-Madjar et al. (1996), 'deducting lines of sight velocity structure for 
D/H evaluations from ionized species could be extremely misleading'. But we fix Za(Silll) = 
0.701024 as a more or less arbitrary reference radial velocity at which Vj = 0. 

The estimated parameters for a few adequate RMC profile fits are listed in Table 1. The 
derived deuterium abundance (D/H) ~ 4.4 x 10~^ is about 4-7 times smaller than the limiting 
values of 1.8 — 3.1 x 10~^ found by Webb et al. (1997b) in the microturbulent model excluding 
the Si III line. To illustrate our results, we show in Figs. 2 and 3 H+D Lya profiles for the two 
calculations with the lowest and the highest D abundances found in the mesoturbulent model 
(D/H = 4.111 X 10~^ and 4.755 x 10~^, respectively). They are shown by the solid curve, whereas 
the open circles give the experimental intensities. The residuals shown in Figs. 2b and 3b by 
filled circles are normally distributed with zero mean and unit variance. 

To control the f (s)-configurations estimated by the RMC procedure, we calculated profiles 
for the higher order Lyman lines and the shape of the Lyman-limit discontinuity and then 
superposed them to the corresponding part of the Q 1718+4807 spectrum. One expects to find 
a displacement between the observed and the calculated steep Lyman edge if the velocity field 
structure is far from the truth. The results are shown in Figs. 2c and 3c where again open 
circles correspond to the observed intensities and the computed spectra are shown by the solid 
curves. We do not find any pronounced discordance of calculated and real spectra. On the other 
hand, the spectral resolution of 2.95 A (FWHM) is not sufficient to follow a fine velocity field 
structure within the Za = 0.701 absorber. 
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Figure 2: Observations and RMC fits for Q 1718+4807. (a) - HST/GHRS data (open circles) 
and calculated profiles for HI (dashed curve), DI (dotted curve) and H+D (solid curve). The 
spectral resolution corresponds to a FWHM = 0.1 A. Shown is an acceptable solution of the 
inverse problem with D/H = 4.111 x 10-^ N(HI) = 1.768 x lO^^^ cm-2, Tkin = 1-51 x 10^ K, 
at = 26 km s~^, and L/l = 3.5 (see Table 1). (b) - residuals e in units of CTnoise (see text), 
(c) - lUE spectrum (open circles) and fit (solid curve) with a spectral resolution corresponding 
to a FWHM = 2.95 A. 
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Figure 3: As Fig. 2 but for D/H = 4.755 x 10"^ N(HI) = 1.759 x lO^'^ cm~^ Tkin = 1-75 x 10^ 
K, at = 40 km s"^ and L/l = 2.8 (see Table 1). 
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Figure 4: Probability density functions p{v) for the velocity components parallel to the line of 
sight for two acceptable RMC solutions shown in Fig. 2a (dotted line histogram) and in Fig. 3a 
(solid line histogram). Both are blue-shifted by about —20 km s"'^ with respect to ^^(Silll) and 
have asymmetric shapes. 



The derived t'(s)-configurations are not unique. Table 1 demonstrates the spread of the 
apparent rms turbulent velocities which range from 18 up to 40 km s~^. It is worthwhile 
to emphasize once more that the projected velocity distribution function (p.d.f.) may differ 
considerably from a Gaussian. Fig. 4 shows an example of such distortions caused by a poor 
statistical sample (i.e. incomplete averaging) of the velocity field distributions for the two cases 
of the lowest and the highest D/H ratios from Table 1 {at = 26 and 40 km s~^, respectively). 
Both p{v) distributions are asymmetric. This is the main reason why the absorption in the blue 
wing of the HI Lya line may be enhanced without any additional HI interloper (s). 

How to test the validity of the RMC results ? 

We show that the H+D measurements from two absorption-line systems do not require either 
a low or a high value of D/H. Both systems at z = 0.701 toward Q 1718-^4807 and at z = 2.504 
toward Q 1009-1-2956 may be modeled within the framework of the generalized radiative transfer 
theory with D/H in the range 4.1 — 4.6 x 10~^. The fitting procedure leads in this case, however, 
to asymmetric velocity field distributions. This fact may be used to test the RMC results by 
additional observations of higher order Lyman lines. 

Indeed, ifp{v) is asymmetric, this will show up in the profile shapes of the higher order Lyman 
lines. For example, for the physical parameters listed in Table 1, the effect becomes visible 
starting from Ly-4 [the Ly-a, -/?, -7 lines are insensitive to the asymmetry of p{v) due to their 
high optical depth]. Fig. 5 shows simulated spectra (convolved with a Gaussian instrumental 



Table 1: Cloud parameters derived from the Lya profile by the RMC method 
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profile of FWHM = 0.1 A) for Ly-4 and Ly-12 using the same p{v) distributions depicted in 
Fig. 4 - dotted curves for at = 26 km s""*^, and solid curves for at = 40 km s~^. As seen, 
line shapes may vary depending on the velocity field structure that makes them asymmetric in 
general. As for microturbulent profiles shown in Fig. 5 by the short dash curves, they remain 
symmetrical. The higher order Lyman lines may be used unambiguously to constrain the velocity 
field structure of the absorption system. Therefore high resolution observations of the whole 
Lyman series become crucial for the D/H measurements. 

Comparison of the primordial abundances of the light elements 



The revealed mean value of (D/H) ~ 4.4 x 10 corresponds to the mass fraction of He 

p ' 



Yp ~ 0.245, the ^Li primordial abundance of ~ 2.1 x 10 ^°, and 77 ~ 4.4 x 10 ^^ 



The estimated Yp and Li/H values lie within the measured regions 0.240 < Yp < 0.246 
(Izotov, Thuan, Lipovetsky 1997) and 1.5 x 10"^° < '^Li/H < 2.0 x 10"^° (Bonifacio k Molaro 
1997). The measurements of Yp may contain, however, some systematic uncertainties. For 
instance, Fields et al. (1996) and Olive et al. (1997) give the value Yp = 0.234 ± 0.003(stat.) ± 
0.005(syst.) and note that deducing abundances from linestrengths under the usual assumptions 
about extragalactic HII regions 'could introduce a significant systematic error', the magnitude 
and the distribution of which are 'not well understood'. 

The above baryon-to-photon ratio yields the current baryon density Vtb ~ 0.03, assuming 
a Hubble constant Hq = 70 km s^^ Mpc^^ and the temperature of the cosmic background 
radiation Tcbr = 2.726 K. This value of Qb being larger than the observed density of baryons 
in stars and hot gas ri^^obs — 0.01 (Fukugita et al. 1996) implies that most baryons are in the 
hot intergalactic gas. 

Conclusion 

It is not unlikely that the same spectra may yield different column density estimations. 
Model dependence of the D/H measurements has been demonstrated within the framework 
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Figure 5: The HI Ly-4 and Ly-12 mesoturbulent spectra for at = 40 km s^^ (solid curves) and 
at = 26 km s^^ (dotted curves). The corresponding physical parameters are listed in Table 1. 
Short dash curves show microturbulent profiles calculated for the mean N(HI) = 1.738 x 10^^ 
cm~^ and the mean Doppler parameter 6(HI) = 25.5 km s~^ of the data by Webb et al. The 
spectra are convolved with a Gaussian instrumental profile of FWHM = 0.1 A. 
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of microturbulent models by Wampler (1996) and on the basis of mesoturbulent models in a 
series of our papers (see LKT and references therein). Macroscopic motions in QSO absorption 
systems have been repeatedly confirmed by direct observations (e.g. Songaila et al. 1995). Our 
analysis shows that such large scale motions may affect the velocity field configuration along 
the line of sight and the relative positions of the DI and HI lines (Levshakov & Takahara 1996) 
making their unambiguous interpretation difficult. We stress again that the reliability of the 
interpretation of the D/H measurements is determined by two factors : (a) improvements in the 
detection equipment and (b) advances in the theory of line formation in turbulent media. 
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